On the Stacking Charge Order in NaY^Os 
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We propose a mechanism for the observed stacking charge order in the quarter-filled ladder com- 
pound NaV^Os. Via a standard mapping of the charge degrees of freedom onto Ising spins we 
explain the stacking order as a result of competition between couplings of the nearest and next- 
nearest planes with the 4-fold degenerate super-antiferroelectric in-plane order. 



There has been a great interest in recent years from 
both theorists and experimentalists in the insulating 
transition-metal compound NaA^Os 0. This material 
provides a rather unique example of a correlated electron 
system, where the interplay of charge and spin degrees of 
freedom results in a phase transition into a phase with co- 
existent spin gap and charge order. NaV20s is the only 
known so far quarter-filled ladder compound Q. Each 
individual rung of the ladder is occupied by single elec- 
tron which is equally distributed between its left/right 
sites in the disordered phase (see Fig. [I}. At T c — 34 K 
this compound undergoes a phase transition when a spin 
gap opens, accompanied by charge ordering 
The experimentally observed two-dimensional (2D) long- 
range charge order in NaA^Os (the afe-plane order 
is shown in Fig. ^) is super-antiferroelectric (SAF), as 
we have pointed out recently 0. The theory of spin- 
SAF transition put forward by us is adequate 
in accounting for, even quantitatively, various aspects 
of the transition in NaV20s. It deals, however, with 
a single (ab) plane, leaving aside the question of charge 
ordering along the third (c) direction. The phase transi- 
tion in NaV205 quadruples the unit cell in c-direction 
(the supercell of the ordered phase is 2a x 2b x 4c), 
and the recent X-ray experiments 0, E| revealed peculiar 
ordering patterns in c-direction (stacking order) of the 
super-antiferroelectrically charge-ordered planes. Here 
we present a model which provides the explanation for 
the observed stacking order in NaV^Os. 

An insulating quarter-filled ladder system, with elec- 
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FIG. 1: NaV2C>5: ladders in the afo-plane. In the disordered 
phase each electron is equally distributed between left/right 
sites on a given rung (left panel), while below T c electrons 
(filled black circles) order as shown in the right panel. 



trons localized on the rungs of the ladder, can be mapped 
on an effective spin-pseudospin model, where the Ising 
pseudospin {T x ) represents left /right positions of the 
charge on a given rung, similar to the standard pseu- 
dospin approach to the order-disorder-type phase transi- 
tions 101. The 2D effective spin-pseudospin Hamiltonian, 
able to describe the spin-SAF transition in NaV2 05, was 
given in . Since the present work is concerned with the 
physics of the charge order, we will discuss here exclu- 
sively the Ising sector of the full spin-pseudospin model. 

Within the spin-pseudospin formalism, the afe-plane of 
coupled ladders can be mapped on the effective lattice 
shown in Fig. [2 which in its turn we smoothly map onto 
a more conventional square lattice. It is then easy to 
identify the charge order in the afr-plane (see Figs.^and 
01 below) as the SAF phase of the 2D nearest neigh- 
bor (nn) and next-nearest neighbor (nnn) Ising model, 
shown in Fig. [31 Since the Ising couplings Jjj = ■/□, Ji, <H 
originate from the Coulomb repulsion, we assume them 
to be antiferro (AF), i.e., Jj > 0. SAF is the ground 
state of the 2D nn and nnn Ising model, if [jj 



J1 + J2 > \ Jc 



and J\ 2 > 



(1) 



The SAF state can be viewed as two superimposed an- 
tiferromagnetically ordered sublattices (circled/squared 
sites shown in the right panel of Fig. |3J, and it is 4- fold 
degenerate, since each of these sublattices can be nipped 
independently. 

Some possible patterns of the stacking charge order 
in NaV205, determined from the X-ray experiments 0], 
are shown in Fig. 31 In terms of the effective Ising model 
this translates into the 3D pseudospin ordering patterns 
depicted in Fig. To explain the mechanism of these 
types of order, let us consider the "minimal" 3D nn and 
nnn Ising model with the Hamiltonian 



k.l 



(2) 



where the bold variables denote lattice vectors. The sum 
includes spins on the nn sites coupled via Jp , J3 and on 
the nnn sites coupled via Ji,2,4 (cf. Figs. 13151) . This 
model reduces to a more familiar 3D ANNNI model 0] 
when the diagonal couplings Ji,2 are absent. We will be 
interested in the case of the plaquette couplings satisfying 
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FIG. 2: a&-plane of coupled ladders mapped onto an effective 
2D lattice. A vertical line and a dot represent a single ladder 
and its rung. The Ising pseudospin (up/down) represents the 
position (left/right) of the electron on a given rung. 
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FIG. 3: Couplings on an elementary plaquette of the 2D nn 
and nnn Ising model (left) and an example of the SAF order 
(right). An plaquette on the square lattice corresponds to the 
encircled region in Fig. [2] 



JU, i.e., when the planes are SAF-ordered in the ground 
state, and a frustrating nnn-interplane coupling J4 > 0. 

Note, that the four patterns A-D of the SAF state 
can be divided into two pairs of the "AF counterparts" 
(A,B) and (C,D) in the sense that in the presence of 
an AF coupling between two SAF-ordered planes, the 
counterparts A and B (or C and D) minimize the energy 
without inter-plane frustrations (cf. Fig. • 

From energy considerations one finds that, if J4 > 
IJ3I/2, the model has a 16-fold degenerate ground state 
composed of the SAF-ordered planes stacked with a pe- 
riod of four lattice spacings in z(c)-direction. We will 
denote this 3D order as SAFx4. (Note that the in-plane 
SAF phase itself can take either 2x1 pattern, or 1 x 2.) 
The number of frustrated inter-plane bonds (J3) in the 
SAFx4 phase is 8 per 4 stacked elementary cubes of the 
lattice (in average 2 J3-bonds per elementary cube) . This 
phase can be realized via four 4-fold degenerate stacking 
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FIG. 4: Possible patterns of the stacking charge or- 
der in NaV2 05 from the X-ray experiments (reproduced 
from Ref.0): (BDAC) (and (DBCA), not shown) [top]; 
(DDCC) (and (BBAA), not shown) [bottom]. Each layer 
(in a&-plane) labeled by A-D corresponds to a particular re- 
alization of the 4-fold degenerate SAF state. 



patterns: 

(AABB), type I 

(ACBD), type II 

(ADBC), type III 

(CCDD), type IV. 

4-fold degeneracy of each of these stacking patterns 
comes from translations along the stacking direction, or, 
in terms of the above notations, from cyclic permuta- 
tions inside the angular brackets. So, the two particular 
realizations of the SAFx4 state shown in the left and 
right panels of Fig. [S] belong to the patterns II and IV, 
respectively. 

If J4 < I J3 1/2 the model can have two possible 4- fold 
degenerate ground states SAFxl or SAFx2, depending 
on the sign of J3. If the nn-interplane coupling is AF J3 > 
0, the model's ground state SAFx2 can be realized via 
two 2-fold degenerate stacking patterns (AB) or (CD). 
For a ferromagnetic nn-interplane coupling J3 < 0, the 4- 
fold degenerate ground state SAFxl is simply one of the 
four possible SAF patterns stacked in c-direction. In each 
of the phases SAFxl or SAFx2 there are 4 frustrated 
inter-plane bonds ( J4) per elementary cube of the lattice. 
Thus, to summarize the ground-state phases: 



J 4 > |^3 1/2 : SAFx4 
J 4 < J3/2 : SAFx2 
Ji < -J3/2 : SAFxl 



(3) 
(4) 
(5) 
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FIG. 5: Two particular realizations of the 16-fold degenerate 
ground state of the 3D nn and nnn Ising model © at J4 > 
IJ3I/2. The depicted Ising orders correspond to the charge 
ordering patterns shown in Fig. 0] 



The model (121! ill provides the explanation for the 
charge ordering in NaV2C>5 found in the most recent X- 
ray studies §|, carried out at ambient pressure deep in 
the ordered phase (at T = 13 K). Those authors found 
that the sum of patterns (BDAC) + (DBCA) (i.e., 
(II + III) type) and/or of (DDCC) + (BBAA) (i.e., 
(IV + 1) type), makes the best fit to the scattering data, 
and not a single pattern. This implies that the actual 
stacking charge order in NaV205 accommodates all those 
(degenerate) patterns with stacking faults 0, E|. 

An analysis of stacking faults between the various (to- 
tally 16) faultless patterns shows that they cost different 
energy. For instance, if J3 < the least energetically 
expensive (the energy is (J4 — | Jj,\/2)N 2 with respect to 
the ground state for N x N x N lattice) are faults of the 
type (II or III) • (I or IV), where • indicates the position 
of the fault. The examples are: 

....ACBD • DCCD. .. (6) 
... ADBC • AABB. .. 

The faults of the type (II or III) • (II or III), (I) • (I), 
(IV) . (IV), like, e.g., 

... ACBD •DACB. .. (7) 
... ACBD .DBCA ... 
... AABB .BAAB.... 

cost twice more energy than those of the type (jHJl . There 
are of course many other possible types of stacking faults 
(with higher energies) which we will not discuss. From 
energy consideration, we can conclude that at low tem- 
perature faults of the type © prevail. This result to- 
gether with the experimental findings @ lead us to the 
conclusion that in NaV2 05 all four possible types of 
ordering (I-IV) occur indiscriminately, with minimal- 
energy stacking faults, e.g., (II) • (I) • (III) • (IV). 



Above we assumed J3 < 0, which might appear odd, 
since all Ising couplings originate from the Coulomb re- 
pulsion, i.e., they are antiferro. In fact such ferromag- 
netic coupling is an effective coupling replacing interac- 
tions between several charges in order to keep the model 
"minimal" 1 In this sense J4 is an effective coupling as 
well, both J3 and J4 should be viewed as phenomenologi- 
cal parameters whose signs and relative strengths are cho- 
sen to agree with experiments. We can however suggest a 
simple mechanism resulting in an effective J3 < 0. Let us 
consider two nn ladders (two chains in terms of the effec- 
tive lattice shown in Fig. 0) separated by a lattice unit c 
in the stacking direction. On a single plaquette in the foe- 
plane, let us take into account the nn couplings J\ (along 
(0, 1, 0)-direction) and J3 (along (0,0, Indirection), with 
(Ji > J3), in addition to a nnn-coupling Jd along pla- 
quette's diagonals (along (0, ±1, Indirections). A simple 
analysis shows that two nearest stacked chains tend to 
align ferromagnetically in c-dircction when Jd > J3/2. 
We can therefore take J3 = J3 — 2 Jd as the single effective 
nn coupling, as shown in Fig. Considering the actual 
lattice distances in NaX^Os and the distance-dependence 
of the Coulomb repulsion one finds J3 < 0. 

Experiments on NaV2C>5 under pressure 13] show that 
the in-plane SAF charge order is robust and does not 
change, while at the pressure P c = 0.92 GPa a transi- 
tion of the ground state charge order from SAFx4 into 
SAFxl occurs. Thus, in NaV2C>5 the pressure depen- 
dence of in-plane couplings is "non-critical", i.e., condi- 
tion Q is satisfied, while the ratio n = J4/J3 of the cou- 
plings between planes is more sensitive to pressure. As 
follows from (|3I5|) this ratio reaches the frustration point 
k c = —1/2 at P c , and then the ground state changes. 2 

Another very interesting feature of the charge ordering 
under pressure in NaA^Os is the existence of a region 
with numerous higher-order commensurate superstruc- 
tures SAFx^ (where m,n are integers), i.e., a devil's 
staircase region above P c on the temperature-pressure 
plane |l3j . Ohwada and co-workers noticed resemblance 
between the experimental phase diagram and that of the 
3D ANNNI model. (For reviews on that model see, e.g., 



1 For single layer a "truly minimal" model, still having the SAF 
phase, would be that with /□ = 0. In this case it consists of two 
decoupled superimposed (AF) Ising lattices. Then our model J2J 
with nn and nnn interactions between layers reduces to two iden- 
tical inter-penetrating decoupled 3D ANNNI models, cf. Figs. 
I Let us call this limit the A(g)A model for brevity. 

2 In connection to what was said before about J3 , we should point 
out that experiments do not rule out J3 > at ambient pressure 
since SAFx4 order is insensitive to the sign of J3, as long as 
is satisfied. Then J3 must decrease under pressure such that 
at P c it is already ferromagnetic and ft = k c . We find such a 
possibility rather exotic, and will not consider it here. Note also 
that if J3 > at ambient pressure, the above analysis of energies 
of stacking faults should be modified accordingly. 
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Refs.mOSl.) However 

an explanation of the charge or- 
dering in NaV2C>5 in the framework of the ANNNI model 
is incorrect, as the latter cannot have the in-plane SAF 
order in principle. The minimal model to reproduce the 
observed experimental results is what we call the A(gA 
model (see footnote^): two superimposed square lattices 
create the observed overall SAF order in the individual 
layers, while the period of the charge order in the stacking 
direction is the same as in a single 3D ANNNI model. 

Considering the geometry of the original NaV20s lat- 
tice, one can see that Ja is indeed rather small in com- 
parison to Ji, 3 so the minimal A<£>A model appears to 
be adequate for the description of the charge ordering in 
that compound. 

The more complicated model J5J) with J n ^ and 
© satisfied, is expected, from mean-field considerations, 
to also show a sequence of commensurate phases origi- 
nated from the frustration point k c = —1/2 separating 
the SAFx4 and SAFxl ground states, similarly to the 
3D ANNNI model. However to substantiate this sugges- 
tion, a separate study of the temperature phase diagram 
of © is warranted. 

So far we have been discussing the charge in NaV20s in 
terms of Ising models. The full spin-pseudospin Hamilto- 
nian of the problem is more involved |jj , since it includes 
also a transverse field in the Ising sector plus coupling of 
the charge (Ising pseudospin) to the spin degrees of free- 
dom. The present 3D extension of the Ising sector can be 
treated along the lines of our earlier analyses 0, 13 , re- 
sulting in simultaneous appearance of the charge (Ising) 
order and spin gap. 

Conclusions: We propose a mechanism for the stack- 
ing charge order in NaA^Cv It is a result of competi- 
tion between couplings of the nearest and next-nearest 
planes with the 4-fold degenerate SAF in-plane order. 
The simplest effective model resulting in the observed 
charge ordering patterns consists of two decoupled inter- 
penetrating 3D ANNNI models (the A® A model). 

3 The same arguments suggest a weak diagonal J2, which is, how- 
ever, reinforced effectively due to the bilinear spin-pseudospin 
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